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Abstract 

The B-factories results provide an impressive confirmation of tlie Standard Model (SM) 
description of flavor and CP violation. Nevertheless, as more data were accumulated, 
deviations in the 2.5-3.5 a range have emerged pointing to the exciting possibility of new 
CP-odd phase(s) and flavor violating parameters in B-decays. Primarily this seems to 
be the case in the time dependent CP asymmetries in penguin dominated modes (e.g. 
B — (f){r)')Ks). We discuss these and other deviations from the SM and, as an illustration 
of possible new physics scenarios, wc examine the role of the Top Two Higgs Doublet 
Model. This is a simple extension of the SM obtained by adding second Higgs doublet 
in which the Yukawa interactions of the two Higgs doublets are assigned in order to 
naturally account for the large top-quark mass. Of course, many other extensions of the 
Standard Model could also account for these experimental deviations. Clearly if one takes 
these deviations seriously then some new particles in the « 300 GeV to « few TeV with 
associated new CP-odd phase(s) are needed. 



1 Introduction 



The spectacular performance of the two asymmetric B-factories is a triumph of accelerator 
science. Both machines appreciably exceeded their designed luminosities and presently have 
delivered over 1 ab~^ of data [1]. 

On the one hand, the first crucial result is a striking confirmation of CKM-paradigm [2] of 
flavor and CP violation. It is clear that the CKM phase provides the dominant explanation for 
the observed CP violation to an accuracy of about 10-15%. This strongly suggests that new 
physics most likely can only show up as a perturbation requiring accurate measurements and 
precise theoretical calculations. 

On the other hand, many B-factory results indicate interesting deviations from the SM. One 
of the most compelling hints of new physics are the measurements of the time-dependent CP 
asymmetries in penguin dominated modes that turned out to be systematically smaller than 
the SM expectation. While the calculation of these asymmetries requires to keep under control 
long distance QCD effects, the QCD factorization as well as several other approaches shows 
that some of the modes are extremely clean {i.e. (f)Ks, Tj'Kg and Kg Kg Kg final states). The 
magnitude of the deviation ranges from 2.5(7 to about 4a depending on how one chooses to 
compare. The amplitudes for these decays arc dominated by penguin {i.e. short distance) 
contributions: hence, deviations in these CP asymmetries are expected and quite natural in a 
very wide class of new physics scenarios. It is therefore extremely important to follow this issue 
very closely. 

Unfortunately a sizable reduction of the experimental errors on these asymmetries requires 
significantly greater statistics and is bound to be slow: the projected doubling of the integrated 
luminosity by the end of 2008 is unlikely to resolve this issue in a decisive fashion. From a the- 
oretical point of view, it would be extremely desirable to reduce the experimental uncertainties 
at the 5% level because a Standard Model irreducible pollution of a few percents is expected. 
The needed luminosity for this important enterprise may have to await the advent of a Super-B 
factory [3-5]. 

In addition to this hint for new physics there are several other measurements that deviate 
sizably from the respective SM expectations. Among those we have the muon anomalous 
magnetic moment, difference in direct CP asymmetries m. B ^ Kir decays and the tension 
between \Vub\ and the time-dependent CP asymmetry in B ^ J/ipRg. In this paper, we 
present an extensive discussion of several important experimental hints for deviations from the 
SM. As an illustration we study how a simple and well motivated extension of the SM (the Top 
Two Higgs Doublet Model [6-8]) can handle the experimental data. 

Needless to say, the deviations seen in B decays [9] and some other aspects of flavor physics, 
may also be accountable by many other extensions of the SM; for example supersymmetry [10], 
a fourth family [11], a Z-pcnguin [12], warped flavor-dynamics [13] etc. Clearly, the key features 
of any beyond the Standard Model scenario that is to account for the experimental deviations 
in B-physics and other flavor physics that are being discussed here are that it has to have at 
least one new CP-odd phase and new particles in the 300 GeV to ^ few TeV range. Much 
more experimental information is required to disentangle the various possibilities. 
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In Sec. [2] we present the list of problematic measurements that we consider and summarize 
them in a pull table. In Sec. [3] we give a short overview of the Top Two Higgs Doublet Model 
(T2HDM). In Sec. |4]we perform a chi-squared analysis of the T2IIDM and show how present 
experimental results, using observables that are relatively clean, constrain its parameter space. 
In Sees. [5] and |6] we present details of the calculation of T2HDM contributions to various 
observables. A brief summary and outlook is given in Sec. [7j 

2 Possible hints for deviations from the SM 

In this section we summarize some of the experimental problems that have surfaced in the past 
few years connected with the Standard Model picture of flavor physics. In particular we focus 
on the tension between the measured time dependent CP asymmetry in i? — »• J/ipRg and the 
rest of the unitarity triangle fit, the discrepancy between CP asymmetries in 6 — > sss {e.g. 
B {(f),'r]')Ks) and b ccs [B J/ipKs) transitions, the difficulties in describing the CP 
asymmetries in the decays B^ K^tx^ and B~ K'n^, the anomalous magnetic moment of 
the muon and the forward-backward asymmetry in Z ^ bb. 

a^K- Standard Model prediction vs direct measurement 

The standard unitarity triangle fit, with the inclusion of the constraints from |\4b/V^6|, Sk-, 
AMb^ and AM^^ predicts a^K = sin(2/3) = 0.78 ± 0.04. Here, we used a simple fit 
in which we use the inputs given in Table [l] and assume that all errors are gaussian (this 
means, for instance, that we combine systematic and statistical errors in quadrature). The 
direct determination of this asymmetry via the "gold - plated" ipK^ modes "yields [14] = 
0.675 ± 0.026 and deviates from the SM prediction by about two standard deviations. In Fig. [T] 
we show the SM fit of the unitarity triangle in the (p, f/) plane and the a^K constraint is 
superimposed. From the figure it is clear that this effect is mainly due to the conflict between 
a^K and \ Vub\- Note also that the former measurement is clean from hadronic uncertainties and 
the latter uses basically a tree-level process. 

In order to test the stability of this 2a effect, it is useful to entertain a scenario in which 
the errors on |\4;,/V^fe| and on the SU(3) breaking ratio obtained by lattice calculations are 
increased. Increasing 5|Ki6/V^6| = 10% (from about 7%) and 6^s = 0.06 (from 3-4%), the 
prediction for sin2j3 does not change much: we flnd sin(2/5) = 0.78 ±0.05. It is also interesting 
to consider the impact of the very recent lattice determination of Bk presented in Ref. [15]: 
using Bk = 0.765 ± 0.017 ± 0.040, the fit gets slightly worse and the prediction for sin 2/5 
reads 0.76 ± 0.035. The conclusion of these exercises is that the strain between the direct 
determination of a^x and the rest of the unitarity triangle fit is quite solid [16]. 

Time— dependent CP asymmetries in 6 — ^ sss modes 

Within the SM, the CP asymmetries in penguin dominated 6 — >■ s transitions such as (pKg and 
rj'Ks are equal to sin(2/5) up to penguin polluting effects, that are expected to be fairly small in 
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Figure 1: Unitarity triangle fit in the SM. The constraints from iKib/Kfel, £k, ^MbJ AMb^^ 
are included in the fit; the region allowed by a^K is superimposed. 

these modes [17-19]. The calculation of matrix elements of penguin operators is an intrinsically 
non-perturbative task, and it has been recently studied using many different approaches [20-23] . 
These studies show that while a precise calculation of hadronic uncertainties is very difficult, 
at least three cases, namely fj'Kg, (pK^ and KgRsKg [24] are notably clean with only a few 
percent contaminations. In many other cases rough estimates (see for instance Refs. [19-21]) 
suggest hadronic uncertainties to be less than 10%. For example, Ref. [21] quotes a^//^ — a^x = 
0.01 ± 0.01 and a^x — o-ipK = 0.02 ± 0.01. The measurements of the time dependent CP 
asymmetries in the rj' and (p modes, yield ar^'K = 0.61 ± 0.07 and a^K = 0.39 ± 0.18. The latter 
deviates from the SM prediction a^ss = 0.78 ± 0.04 at the two sigma level. 

It is rather curious that all the time dependent CP asymmetries in 6 — > sss have been 
measured to be somewhat smaller than the B J/ipKg asymmetry. If we naively compute 
the average of the CP asymmetries in all the h — sss modes, even though only three of the 
modes are very clean and others may have 0(10%) uncertainties, one then finds [asss]average = 
0.52 ± 0.05 with a deviation of about 4cr from the SM prediction and about 3(j from the value 
directly measured by the ipK^ method. 

For the sake of completion, we also note that just averaging over the three clean modes 
gives ttciean = 0.57 ± 0.06. Since in so far as the SM is concerned, sin 2/3 may be measured 
either by these three clean penguin modes or by the J/tpKs modes, the best "SM" direct 
measurement of sin 2/? is given by the average over the {J/ilj,(f),7]' , KsKs)Ks modes: we thus 
find sin 2/5 = 0.66 ± 0.02 which is again about 2.5cr from the SM prediction of 0.78 ± 0.04. 

CP asymmetries in B — ^ Ktt 

The QCD-factorization predictions for the individual CP asymmetries in i? ^ Kit decays [30, 
31] are extremely sensitive to non-factorizable (hence model dependent) effects and cannot be 
used to directly constrain the SM. Luckily it turns out that, in the calculation of the difference 
between the CP asymmetries in K^n^ and B^ K~tt^, most model dependent 

uncertainties cancel and the QCD-factorization prediction is quite reliable. The magnitude 



3 



iKb/Kfel = 0.1036 ± 0.0074 [25] 


exp _ 
— 


(2.280 ±0.013) 10-3 


Am2;P = (17.77 ± 0.10 ± 0.07)ps-i [26] 


cxp 


= 0.675 ± 0.026 


Am^7 = (0.507 ± 0.005)ps-i 


Bk = 


0.79 ±0.04 ±0.08 [27,28] 


= 1.210^°:°^^ [29] 







Table 1: Inputs that we use in the unitarity triangle fit. 



of this cancelation is apparent in the comparison between the predictions for the individual 
asymmetries and for their difference. The results of Ref. [31] read: 

Acp{B-^K-^') = (7.lll:™l°n?) % (1) 
AcpiB'^^K-n^) = (4.5l}:}l|^l°|l«:f) % , (2) 

where the first error corresponds to uncertainties on the CKM parameters and the other three 
correspond to variation of various hadronic parameters; in particular, the fourth one corre- 
sponds to the unknown power corrections. The main point is that the uncertainties in the two 
asymmetries are highly correlated. This fact is reflected in the prediction for their difference; 
we find: 

AAcp = Acp{B- K-TT^) - Acp{B^ R-n^) = (2.5 ± 1.5)% . (3) 

In evaluating the theory error for this case, we followed the analysis presented in Ref. [31] and 
even allowed for some extreme scenarios (labeled S1-S4 in Ref. [31]) in which several inputs 
are simultaneously pushed to the border of their allowed ranges. The comparison of the SM 
prediction in Eq. ^ to the experimental determination of the same quantity [14] 

AA^^I = (14.4 ±2.9)% , (4) 

yields a 3.5a effect. 

Muon g — 2 

The muon anomalous magnetic moment has been thoughtfully investigated in the literature. 
The most up-to-date calculation of the SM prediction suffers from model dependent uncertain- 
ties in the calculation of the light-by-light hadronic contribution; nevertheless, all the estimates 
(see, for instance, Ref. [32] for a collection of results) point to a SM prediction that is lower 
than the experimental measurement by about three sigmas. The inconsistency between the 
extraction of the hadronic contribution to the vacuum polarization from r and e~^e~ data is 
still an open question. We note, however, that the use of the former requires model dependent 
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assumptions on the size of isospin breaking effects; for this reason, most analyses prefer not to 
include r decay data. The most recent theory estimate is [33] 

a™ = 116591785(61) x 10"^^ , (5) 
while the present measurement is [34,35]: 

af^ = 116592080(63) x 10"^^ . (6) 
The discrepancy is at the 3a level. 

Forward— backward asymmetry in Z — ^ bb 

The LEP measurement of the forward-backward asymmetry in Z ^ bb reads ^4°'^^ = 0.0992 ± 
0.0016. The discrepancy between this experimental result and the central value of the SM fit, 
{A^jr'l^)sM = 0.1038 is at the 3 sigma level. Care has to be taken in interpreting this result 
because the indirect determination of A^jj^ from the forward-backward Left-Right asymmetry 
(Af,) is compatible with the SM prediction at 1 sigma. 

Overview: the pull table 

Let us give a global view of the status of the Standard Model by collecting most measurements 
sensitive to the flavor sector and their deviation from the corresponding SM prediction^ Note 
that several of the entries indicate deviation from the SM in the 2.5 - 3.5 cr range. 



Observable 


Experiment 


SM 


Pull 


B{B X,7) X 10^ 


3.55 ±0.26 


2.98 ±0.26 


±1.6 


B{B TUr) X 10^ 


1.31 ±0.48 


0.85 ±0.13 


±0.9 




17.77 ±0.12 


18.6 ±2.3 


-0.4 




0.675 ±0.026 


0.78 ±0.04 


-2.0 




0.39 ±0.18 


0.80 ±0.04 


-2.2 




0.61 ±0.07 


0.79 ±0.04 


-2.0 




0.51 ±0.21 


0.80 ±0.04 


-1.3 


0'{<i)K+-q'K+KKK) 


0.57 ±0.06 




-2.9 


^{(PK+Ti'K+KKK+ipK) 


0.66 ±0.02 




-2.6 



*See Sec. 5 for a detailed discussion of the various observables 
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L £>t)o J llcLiVcclVcl ctgfc! 


0.52 ± 0.05 




-3.7 


Ar,/r, 


0.27 ±0.08 


0.147 ±0.060 


±1.2 


AAcp 


0.144 ±0.029 


0.025 ±0.015 


±3.6 




1.16592080(63) 


1.16591785(61) 


±3.4 


.0,6 
^/6 


0.0992 ± 0.0016 


0.1038 


-2.9 


\Vub\ X 10^ 


4.31 ±0.30 


3.44 ±0.16 


±2.6 



In the next section we will introduce a particular new physics model, the two Higgs doublet 
model for the top quark (T2HDM), and see how well it can accommodate the above mentioned 
deviations. 



3 The two Higgs doublet model for the top quark 

The T2HDM is a special case of type-Ill 2HDM. It was first proposed in Ref. [6] and subse- 
quently analyzed in Refs. [7,8,36]. In this model, one of the Higgses has only interactions in- 
volving the right-hand top, while the other one couples to the remaining right-handed fermions 
(but not to the top). The main motivation for this model is to give the top quark a unique 
status, thus explaining in a natural way its large mass; hence large values of tan jdn (the ratio of 
the vev's of the two Higgs fields) are preferred. As we will see in the following, a consequence of 
the pecuhar structure of the T2HDM is that the model contains two additional flavor changing 
complex couphngs on top of the standard 2HDM parameters. 
The Yukawa interactions of the quarks with the Higgs flelds are: 

Cy = -Qii/iF,ciR-gi^iy,i(i2W-QL^2>;i^'W±h.c. , (7) 

where are the two doublets, = ia"^!!*, Y^^d are Yukawa matrices, 1*^^^-' = diag(l,l,0), 
and l*-"'^ = diag(0,0, 1). After the clcctrowcak symmetry breaking, the neutral components of 
Higgs doublets receive two independent complex vev's, vi/ \/2 = ve'''^^ cos (3h/V^ and V2/ a/2 = 
ve^'^^ smPn/V^, whose ratio is tanPn = \v2/vi\. 

The quark mass matrices in the mass eigenstate basis are: 

m, = Dli^^Y^jD^, (8) 
mu = t/i ± C/« , (9) 
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where mu,d are diagonal, Ul,r and Dl^r are unitary matrices and V — Uj^Dj^ is the CKM 
matrix. The charged and neutral Higgses interactions read: 



J^Y — 



-ULVrriddR / - ur m^V / + EV 



di + h.c. 



^2 



u 



VrUdPR + m^VPL) (G+ - tan (3hH+) + EVPL(tan (3h + cot (3h)H^ 



MPlrH^ + PLRG^)dR + ur{PSlH^ + PiLG^)dL + h.c. 



d + h.c. 
(10) 



jCy = -dLTUdd 



R- 



0* / H^* 
1= — ul rriu ^—F= + E 



rrO* 

-"2 



0* 



Ur + h.c. 



6/2 tan /3i/ 



2m 



w 



(- , _ \ /i° sin Qiij - cos ctij .(- , _ \ .(, 

I dL-niddR + ulitIuUr ) —z h ^ I dLTUd^^dR + ulttIu^sUr ) 



E'Wi? r^— «(l + COt PHjUL^^-fQURA" 

sm 



+ h.c. 



^ d^Pf + + t7,PfA")dR + udP^^'h' + P„-"if" + ^75Pr ^")«« + h.c. , (11) 

where E = tjIuUrI^'^^Ur. The would be Goldstonc boson G^, the charged Higgs the heavy 
and light scalars and and the pseudoscalar are given by: 



1 
P 



GO 



1^2! 



with 



cos a; — sm a; 
sin Lo cos a; 



The explicit expressions for the charged Higgs couplings are: 



^LR — 



pH 
^RL 



92 



tan Ph V rud , 



LR 



RL 



tan 13h f(l + tan'2 p^)T) - m^V , 

92 



\/2mw 
92 



V rud , 



V2m 



ruu V . 



(12) 
(13) 

(14) 

(15) 
(16) 
(17) 

(18) 



w 
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The explicit expressions for the neutral Higgs couplings are: 

^0 g2 tan (3h 



P: 



2m 



m 



sm an 
w \ sm (Jh 
g2 tan I3h ( cos an 



^cos{aH 



H 



AO 



P: 



2mw 
§2 tan jSu 

g2 tan Ph 
2mw 



mu 



sm ph 
1 + tan^ (3h , 
tan^ /3_ 



sin^ J 

sm{aH - Ph) 
sin^ Ph 



H 



rrid 



sman 



sin P^ 



H 



g2 tan Ph cos an 

md . 

sm Ph 



P. 



AO 



2m\Y 
§2 tan Ph 



2m 



-rrid 



w 



(19) 
(20) 
(21) 
(22) 
(23) 
(24) 



From the definition of S it is clear that only the third row of the matrix Ur is relevant up to 
an overall phase (i.e. we can take {Ur)^^ real). Taking into account the unitarity constraint, 
it follows that S depends on only 4 real parameters. This statement can be explicitly verified 
by employing the most general parametrization of a unitary matrix: U = P1VP2, where Pi are 
diagonal phase matrices and is a unitary matrix that depends on three angles and a single 
phase (e.g. it is CKM-like). The third row of this matrix can always be written as: 



R 



* 



* 

* 



(25) 



r/|2 



where ^ and ^' are complex parameters with |^*^')| < 1. In models based on dynamical top- 
condensation [37,38] and top-color [39,40] the parameters are naturally of order = m^/mt 
(see also Ref. [41]); for this reason we introduce new parameters ^^'^ = ectt^'^ with ^^'^ = 0(1). 
Neglecting terms proportional to the u-quark mass, the matrix S reads: 



E 
mt 



I 















2 C* 



\^ct^^^-\m^-\i?) e.,^Vi-ieiVi-il? (1 
/ 



leivi-iep 
iep)(i-ief) 



eli^ 
Vct^ ^cti 1 



x+0 e 



mt 



(26) 



From Eq. (10) we find the following charged Higgs interactions between right-handed up quarks 



and left-handed down quarks: 



(72 mc tan Ph 
\/2mw 



e Vtd e Vt, 

CVtd-V,d r^ts-Ks CVtt-V,h I . (27) 

V Vtd cot^ pn/ect + ectiVcd Vts cot^ Pnl^ct + erf^K^ Vth cot^ Pnldct , 
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Figure 2: Dependence of Xmin A ''t^h±, tan/Jj^, |^|, ip^, sin 2/?, a and 7. For each value of 
the parameter on the x-axis, we minimize the chi-square with respect to all the others (including 
p and f]. 



In particular t^qL (q=d,s) interactions are dominated by the ^ term for tan Ph > 10. In 
conclusion, the parameters of the models are: tenaPff, an, ^H±y ^ho, ^h": ^% C 

Finally let us comment on the renormalization scheme of the quark masses that appear in 
Eqs. (10) and (11). In the calculation of the additional matching conditions to various operators 
we integrate out the charged and neutral higgses at some high scale /iq ~ 0{mw, fnt); therefore, 
it is most natural to evaluate all the relevant couplings in the MS scheme at the high scale. 
This observation has a very strong impact on the phenomenology of the T2HDM because of 
the strong renormalization scale dependence of the charm quark: m^^ ^ (po) / m^°^'^ ~ 0.45. 



4 Global analysis 

In this section we present the results of global fit of the T2HDM parameter space. Here 
we just focus on the outcome of the fit and investigate how well the T2HDM can answer to 
the problems we collected in Sec. [2| A detailed discussion of the various observables that we 
consider is given in Sees. [5]and|6|, in which we collect the experimental data and the analytic 
formulae required to calculate T2IIDM effects. In those sections we also show the separate 
impact of each observable on the T2HDM parameter space. 

We classify the various observables we consider according to whether neutral Higgs exchange 
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Figure 3: Contour plots corresponding to Xmin < 1 in the (m//±, tan/5//) and v^g) planes. 
For each point on the contour, we minimize with respect to all other variables. The dashed 
and dotted contours correspond to = (1,2) and tan/5j/ = (30,50) for the left and right plot, 
respectively. 

contributions are relevant or not. In the latter case, the parameter count of the model is reduced 
to the sole tan/3/f-, m//±, ^ and Observables insensitive to the neutral Higgs sector of the 
T2HDM include: rare B decays {b 57, b —>■ si^i~ , B ru), neutral meson mixing {K, Bd, 
Bg, D), various CP asymmetries (time-dependent asymmetries in 6 ^ ccs and b sss decays, 
asymmetries in flavor specific B decays, direct asymmetries in the B Ktc system) and the 
neutron electric dipole moment (EDM). Among those observables that display some sensitivity 
to the neutral Higgs sector we consider the muon anomalous magnetic moment, Ap and the 
Z ^ bb vertex {Rb and the forward-backward asymmetry Af,). 

In the analysis we focus on the first set of observables and treat separately the ^' = 
and ^' 7^ cases. In fact, the parameter ^' is related exclusively to transitions between the 
first and third generations and impacts only B tu, DD mixing and the neutron EDM, 
while being completely negligible in all other observables. The T2HDM phenomenology of 
observables dominated by neutral Higgs exchanges is very similar to the one of a regular Two 
Higgs Doublet Model and we will briefly summarize it in Sec. |6} 

Our general strategy is to include directly into the fit only processes for which the theory 
error is reasonably under control; once a region of the T2HDM parameter space has been singled 
out, we look at the other observables. 

r = 

As a first step we set ^' = 0. The that we consider includes the following quantities: 
iKift/Kfel, AMbJAMb^, CLtpK, ^K, B Xs7, B — > TP. The resulting function depends on the 
CKM parameters p and 77, and on the T2HDM parameters mH±, tan/3//, ^ = l^le*"^^. 

Note that, without the inclusion of T2HDM contributions, the overall x^^fit in the SM is 
relatively poor (Xmm ~ 6). Once T2HDM effects are included, the fit improves drastically and 
we find Xmin ~ 0- This implies that this set of measurements singles out a clear sector of the 
parameter space not compatible with the T2HDM decoupling limit. In Fig. [T] we show the 
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Figure 4: Contour plot corresponding to Xmin < 1 in the {^',ip^') plane. The rest of the 
parameters are chosen so to minimize the for ^' = 0. In the plot on the right the light 
gray, dark gray and black regions correspond to a neutron-EDM (given in units of 10~^^ e cm) 
smaller than 3, between 3 and 6.3, and bigger than 6.3, respectively. 



unitarity triangle fit in the Standard Model; note, in particular, the tension between the black 
contour and the constraint from a^^ (not included in the fit). In Fig. [2| we show the actual 
dependence of the full on the CKM angles and the four T2IIDM parameters. The 68% C.L. 
intervals that we find are: 



(660^1°) GeV , (28) 
tan/3 = 28+f , (29) 
C > 0.5 , (30) 

n = {imiy, (31) 

p = 0.19 ±0.035, (32) 
f] = 0.38 ±0.03. (33) 

The corresponding ranges for the three UT angles are: 

sin(2/3) = 0.77 ±0.04, (34) 
a = (89 ± 6)° , (35) 
7 = (64 ± 5)° . (36) 

In Fig. [3] we show the correlation between these parameters; the shaded areas correspond to 
Xmin — 1 their projections on the axes yield the corresponding la regions. 

We are now in the position of evaluating how well the T2HDM does with respect to the pull 
table we introduced in Sec. |2l From the outcome of the fit it is clear that the T2HDM can 
easily accommodate the deviations in B ^ X^'j, B tu, a^x and \Vub\. Unfortunately, for 
> 400 GeV, it seems quite difficult to accommodate the effect required to reconcile the 



CP asymmetries in B ^ (rj', (j))Ks with experimental data (see Fig. 13 in Sec. 5.4). The impact 
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of the T2HDM on AMb^ is not very large and is perfectly compatible with the present exper- 
imental determination. Finally we do not find any large contribution to the CP asymmetries 
in B ^ Ktt, hence within the T2HDM the 3.6a observed deviation remains unexplained. 

Let us note that the solution of the \Vub\ vs a^x puzzle is achieved via sizable and highly 
correlated contributions to both a^K and Ek'- after the inclusion of the constraints from B 
Xs'j and B tu, the solution of this puzzle was a crucial bonus that we could not enforce 
(due to the extremely reduced number of parameters that we are considering). 

Finally we point out that T2HDM effects on a^x are caused by large complex contributions 
to the amplitude A{B J/ipKs) and not to the B — B mixing matrix element (i.e. M{^2)- 
Since the former is dominated by the tree-level transition b ccs, any other process controlled 
by this quark-level decay will display similar large effects. This is particularly true for time 
dependent CP asymmetries in Bg decays. The Bs J/ip v' mode, for instance, is based on 
the b — > CCS amplitude, hence, in the naive factorization limit, the T2HDM contributions to its 
time dependent CP asymmetry must be identical to the corresponding ones in i? — J/ipRg. 
Therefore, the above analysis predicts the T2HDM contribution to this asymmetry to be in 
the +10% range. Given that the SM expectation for this quantity is extremely small (the phase 
of the SM Bg — Bs amplitude is about one degree), the measurement of a large enhancement 
in the B —>■ J/ ipK asymmetry is a clear indication for a resolution of the apsiK puzzle via new 
physics in the amplitudes (as it is the case in the T2HDM). 



In order to study the effects of non vanishing we fix the other parameters to the values that 
minimize the just studied; then we include contributions from B tu, DD mixing and 



the neutron EDM (for the latter two, we impose upper limits - see Sees. 5.3 and 5.7 for details). 
Note that without the inclusion of B —>■ tu, the fit for C = ^ favors values of smaller than 1 
(the actual value that we use in the ^' 7^ fit is ,^ ~ 0.8). 

In Fig. |4| we plot the region of {^',(p^') plane for which this new is smaller than 1. 
The main constraint comes from B ru, whose branching ratio is proportional to C,'"^. It is 
interesting to dissect contributions to the neutron EDM: in the right plot in Fig. |4] the regions 
with increasing darkness correspond to a neutron-EDM (in units of 10~^^ e cm) smaller than 
3, between 3 and 6.3, and bigger than 6.3, respectively. 



5 Observables: the charged Higgs sector 

5.1 B ^Xs-f and B XJ+i- 

The experimental world average from the CLEO [42], Belle [43,44] and BaBar [45,46] collabo- 
rations is given by [47]: 

BR{B ^ X,7)s,>i.6GeV = (3.55 ± 0.24^°;?^ ± 0.03) x 10"^ . (37) 
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Figure 5: Plot a. mH± dependence of the branching ratio B —>■ Xs'~f in units of 10"*^. Sohd, 
dashed, dotted and dotted-dashed hues correspond to (tan/3//,^) = (10,0), (50,0), (50,1) and 
(50,-1), respectively. There is no appreciable dependence on The two horizontal dashed 
lines are the experimental 68%C.L. allowed region. The blue region represents the theory 
uncertainty associated to the solid line (similar bands can be drown for the other cases). Plot 
b. Portion of the (tan/3j:^, m^^i) plane excluded at 68%C.L. by the B Xg'j measurement. 
The shaded area corresponds to ^ = 0. The dotted and dashed lines show how this region 
changes for ^ = 1 and —1, respectively. 

The B — i> Xgi^i^ branching ratio has been recently measured by both Belle [48] and BaBar [49] ; 
in the low dilepton invariant mass region, 1 GeV^ < m^^ < 6 GeV^, the experimental results 
read 



B{B ^ X, 
B{B ^ Xs 



fl.493 ± 0.504+°:^^}) X 10-*^ (Belle) 

^-6 



r) = (1.8 ± 0.7 ± 0.5) X 10"^ (BaBar). 
This leads to a world average 

B{B^XJ+e-) = (1.60 ± 0.51) X 10-^ 
The effective Hamiltonian responsible for the transitions b —>■ s'j and b — > si~^i~ is [50] 



n 



■eS 



10 



li=l 



i=3 



where the most relevant operators are 



P7 
Ps 
P9 

Pio 



167r2 
9 



mb{sL(y^''bR)F^y , 



167r 



(38) 
(39) 



(40) 
(41) 

(42) 

(43) 
(44) 

(45) 
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tan (3) 

Figure 6: Plot a. mH± dependence of the branching ratio B XgfJ^fi in units of 10~^. Sohd, 
dashed and dotted hues correspond to (tan/3H,0 = (10,0), (50,1) and (50,-1), respectively. 
There is no appreciable dependence on The two horizontal dashed lines are the experimental 
68%C.L. allowed region. The blue region represents the theory uncertainty associated to the 
solid line (similar bands can be drown for the other cases). Plot b. Portion of the (tan Ph, ^h±) 
plane excluded at 68%C.L. by the B — > Xs/i/i measurement. The shaded area corresponds to 
^ = 0. The dotted and dashed lines show how this region changes for ^ = 1 and —1, respectively. 



The leading order charged Higgs contributions in the T2HDM to the Wilson coefficients 
^7,8,9,10 have been in explicitly calculated in Refs. [7,8,36,51] (see Eqs. (7-15) of Ref . [36]). The 
formula for the new physics contribution to C-j is: 



-1 



i=u,c,t K u a 

2 n D/'„, N I t* +-„„2 



B{yt) + tan^ /3h B{y,) + C tan^ (3h 



mf 



■l^Myc)-el,^B{y,) 



(46) 



where both quantities in square brackets are positive for any choice of tan and mH± , Va = 
m'^/m'jj± and the loop-functions A and B are given in Ref. [36]. 

A numerical formula for the calculation of the B Xg'y branching ratio is given in Ref. [52, 
53]: 



i3(i? -.X,7)t,>1.6GeV 



10" 



2.98 + 4.743 \6Cj\^ + 0.789 I^Cgp + Re( 



-7.184 + 0.612 i) SCr + (-2.225 - 0.557 i) 5Cs + (2.454 - 0.884 i) SCg 5C; 



(47) 



where the leading Wilson coefficients at the high scale are given by Cf°^(/io) = C'i5M(/^o) + 
6Ci and the next-to-leading matching conditions are assumed not to receive any new physics 
contribution, C,-^''(/io) = cI^smIi^o)- The formula above has been obtain by observing that 
using the same numerical inputs of Ref. [54,55] and taking {fic, fib, muo) = (1.5,2.5, 120) GeV, 



the NLO central value of the branching ratio coincides with the NNLO one. Eq. (47) also 
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include an estimate of the new class of power corrections identified in Ref. [56] and of the 
analysis of the photon energy spectrum presented in Ref. [57]. The analyses in Refs. [54,57] 
yield B{B X^'y) = (2.98 ± 0.26) x 10"'*; we will therefore assign a theoretical error of 8.7% 



to the central values calculated in Eq. (47). 



The Standard Model matching conditions and numerical formulae for the calculation of the 
integrated B —>■ Xsi^i~ branching ratios is given in Ref. [50] : 



2.1913 - 0.001655 J(i?io) + 0.0005 I{RwRl) + 0.0535 I{Rr) + 0.02266 liRjRl) 
+0.00496 liRrRl) + 0.00513 X(i?8) + 0.0261 liRgRl) - 0.0118 I{Rg) 
-0.5426 7^(i^lo) + 0.0281 7^(i^7) + 0.0153 n{R7Rlo) + 0.06859 TZiRyRl) 
-0.8554 niRjRl) - 0.00866 n{Rs) + 0.00185 n{RsRlo) - 0-0981 TZiRsRl) 
+2.7008 n{Rg) - 0.10705 TZiRgRlo) + 10.7687 |i?ioP + 0.2889 ji^rT 



+0.00381 |i?8p + 1.4892 li^gl^ 



X 10" 



(4J 



'1.59 + 0.11)10-'' and 



where /2i = Q(/io)/C™(/io). The SM prediction is 5/?(5 
we will assign a theoretical error of 6.9% to the central values calculated in Eq. (48). 

The impact that the B —>■ Xsj and B — > Xgii measurements have on the T2HDM parameter 
space is shown in Figs. [5] and [6j In Fig. [5^ we plot the B — > Xsj branching ratio as a function 
of the charged Higgs mass for various choices of tan Ph and ^. The tan dependence of the 
charged Higgs contributions to Cy is not very strong as it follows from the proximity of the 
solid and dashed curves. The ^ dependence is, on the other hand, much stronger; here we plot 
results for ^ = (1, —1) (other choices of the phase yield in between curves). This can be seen 
explicitly in Fig. |5}d, where we plot the allowed region at 68% CL. in the (tan/3/f, m/^i) plane 
for various choices of ^. Comparison of the plots in Figs. [5] and |6] shows that B Xgii does 
not provide additional constraints on the parameter space. 



5.2 Neutral mesons mixing 

The off-diagonal element of the neutral K-mesons mass matrix is M*^ = {K |Kff| K°)/i2mK), 
where the effective Hamiltonian is 

Kff = ^{^{VZV,,fY.Ca{f^)Qa, (49) 

with 

= i^RdL) {sLdR) (50) 
gf ^ = {sRa^M (sRa^'dL) • 
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mw = 80.426 GeV 


A = 0.818 ±0.012 [27] 


rriK = 0.497648 GeV 


p = 0.197 ±0.031 [27] 


mc(mj = (1.224 ± 0.017 ± 0.054) GeV [58] 


f] = 0.351 ± 0.020 [27] 


mt,poie = (171.4 ±2.1) GeV [59] 


af^{mz) = 0.1182 ± 0.0027 [60] 


sin^ Ow = 0.2312 


mz = 91.1876 GeV 


ml^ = (4.68 ± 0.03) GeV [61] 





Table 2: Numerical inputs that we use in the phenomenological analysis. Unless explicitly 
specified, they are taken from the PDG [25]. 




200 400 600 800 1000 lo 20 30 40 50 

Itlu*^ tan (3) 



Figure 7: Plot a. mH± dependence of the T2HDM contributions to Am^ in units of 10^^ ps^^. 
Solid, dotted and dashed lines correspond to tan/5// = 10, 25 and 50, respectively. There is 
no appreciable dependence on ^ and The horizontal dashed line corresponds to Arri^ < 
0.3 Am^''. Plot b. Portion of the (tan/5//, m//±) plane excluded by the ArnJ^ < 0.3 Am^^ 
constraint. 
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r = 0.985 [62] 


- 1-32 (41)) 


± 0.32 [63] 


Bk = 0.79 ±0.04 ±0.08 [27] 


Jk = 0.159 GeV [27] 


ri2 = 0.57 ±0.01 [64] 




ri3 = 0.47 ± 0.05 [64] 


pLR = -36.1 [65] 


pLR = 59.3 [65] 




pSLL = _i8.i [65] 


pSLL ^ _32_2 [65] 


Am^^P = (5.301 10-3 


) ps-^ 


e^^P = (2.280 ±0.013) 10^=^ 



Table 3: Inputs that we use in the phenomenological analysis of K — K mixing. 



/^^y^ = (0.281 ± 0.021) GeV [66] 
/b,//b, = 1.20 ±0.03 [29] 
mf^(2 GeV) = (0.076 ± 0.08) GeV [67] 
/b, = (0.216 ± 0.022) GeV [29] 

= fsJfBjBjB, = 1.210tni [29] 
Am^J = (0.507 ± 0.005)ps-i 
Am^'^P = (17.77 ± 0.10 ± 0.07)ps-i 


= 5.36675 GeV 
Pl^§^ = -0.89 

P2% = 1-13 
pf;^^ = -0.46 

Pi,h^ = -0.90 

ri^ = 0.55 [62] 


ms, = 5.2794 GeV 
pLji^ = -0.98 

P^,l = 1-24 
pSLL ^ _o.51 
pSLL ^ _Q 98 

a^^P^ = 0.675 ± 0.026 


/d = 0.165 GeV [68] 


mo = 1.8645 GeV 


Bd = 0.78 ±0.01 [69,70] 



Table 4: Inputs that we use in the phenomenological analysis of Bg — Bq and D — D mixing. 



The additional operators Q^^^, Qf^^ and are obtained from Q^^^, Qf^^ and Qf"^ by 

replacing L with R. The effective Hamiltonians that describe B and Bg mixing are obtained via 
the replacements (s, (i) — *■ {b,d) and (s, d) — > (6, s), respectively. The D mixing Hamiltonian 
requires (s,rf) ^ (c,m) and V*gV^^ V^V^f,). 
In the SM only the coefficient C^ll 

receives sizable contributions (in the D meson sector the 
GIM cancelation is more effective due to the smallness of the b quark with respect to the top 
one). 

In the T2HDM there are no tree-level flavor changing neutral Higgs currents involving down 
quarks; hence the Wilson coefficients for K, B and Bg mixing receive non standard contributions 
only through charged Higgs box diagrams. The latter can be found, for instance, in Eq. (A. 11) 
of Ref. [71] g 



The situation is different for what concerns D — D mixing. In fact, from Eqs. (11) and 



(26), it follows that the ulCrS^ {S = h, H, A) coupling is non-vanishing (albeit quite small); 



^We defined the couplings P^j^ ^'i- (^^l in complete analogy to Ref. [71] 
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(a) (b) 




Figure 8: Plot a. mH± dependence of the T2HDM contributions to Ek in units of 10"'^ 
{e^ = epHDM _ ^SM^_ Solid, dotted and dashed lines correspond to |^| = 1 and tanPn = 10, 
20 and 40, respectively. Curves with positive and negative correspond to ,^ = (1,-1), 
respectively. There is no appreciable dependence on The meaning of the blue region is 
explained in the text. Plot b. Portion of the (tan Ph, ^h±) plane excluded by Ek- The shaded 
area corresponds to ^ = 1. The dashed line show how this region changes for ^ = e*'^/^. Other 
choices of the phase yield in-between lines. 



therefore, it induces a tree level contribution to the Wilson coefficient Cf . The charged 
Higgs box diagram contributions are obtained from Eq. (A. 11) of Ref. [71] with the following 
replacements: d ^ u, ^ (Pf )t (for A = LR, RL and B = G, H), V V\ LR ^ RL and 
{ji) — i> (21). Neutral Higgs box diagrams involve the small ulCrS^ coupling and are suppressed 
with respect to the tree level contributions. 



5.2.1 KK mixing 



The K — K mass difference and the measure of indirect CP violation in the K system are given 
by (see for instance Ref. [64]) 



Sk = 



2 Re(Mi^) , 

A{Kl (7r7r)/=o) ^ exp(i7r/4) 
A{Ks (7r7r)/=o) V2AmK 



Im(M/^) 



The expression for M^^ in presence of arbitrary new physics contributions is [62,64,65]: 



?K 



Or p 

12^ 



fj^BKmRm^, A* rjiSoixc) + X*t r]2Ftt + 2XlXlri3So{xc,Xt 



4r 



r~<VLL 
'-'new,K 



, r^VRR , nLRr^LR , pLRr^LR 



/^SLL 1 /^SRR I pS 



K 



riSLL I r^SRR 
^2,K "T ^2,K 



(51) 
(52) 



(53) 



(54) 



I pSLL 

where Aj = V^^^^, Xt = /rri^r, Xc = Ml/my^, the functions 5*0 are given for instance in 
Ref. [64], Tji and r are the QCD correction to Sq^Xi) in the SM, fx is the kaon decay constant. 
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Bk and Pf" = Pf^/ (4772-8^) are lattice QCD determinations of the matrix elements of the 



operators in Eq. (50) [62,65]. The numerical inputs the we use are summarized in Tables |2]|3} 
The KK mass difference receives additional long distance contributions; in the numerical 
analysis we assume that such non-perturbative effects do not contribute to more than 30% of 
the observed mass splitting (i.e. {^mK)NP < 0.00159 ps~^). See Ref. [72] for an estimation of 
these long distance effects in the large limit. An approximate expression for ArriK is the 
following: 



ArriK 



67r 



' 2 , ^2 ["^c("^t)] tan^/^H^ 

r/i mj + — 

^ r 4 mj^± 



(55) 



Imposing the (Ami<-)7vp < 0.00159 ps^^ constraint, we obtain: m//± > 89 (tan/?j^/25)^ GeV. 

The exact numerical impact of the upper limit on ( Ami<-) tvp can be seen in Fig. [7| Compar- 
ison with Fig. [5]d shows that for ^ > 0, this constraint is complementary to -B ^ ^sl- 

The impact of the Ek measurement is shown in Fig. [s] Here we require ^puDM 
Ek range extracted from the standard unitarity triangle analysis. A more correct approach is to 
fit the unitarity triangle in the T2HDM and check whether each given point in the parameter 
space gives an acceptable chi-square. This analysis is presented in Sec. |4j 

We find that the inclusion of the Ek constraint has a very strong impact. Note that, in this 
case, the effect is proportional to ^; hence, B — > Xs7 is still required in the ^ ~ limit. 



5.2.2 BqBq mixing 

The Bq — Bq mass difference is given by [71] 



G%m 



w 



ms^VEflBs, 



(56) 



where F^^'' 



is given by Eq. (541) with the replacement K —>■ Bq. We recalculated the quantities 



Pf" = Pfl {ArjBBBq) using the formulae presented in Ref. [65] and the lattice results of Ref. [73, 
74] . The numerical inputs that we use are collected in Table |4j 

The SM prediction for Am^^ does not depend on the extraction of the CKM parameters p 
and 1]] using the inputs summarized in Table |4| we obtain Am;|^^ = (20.5 ± 3.1) ps~^. Note 
that, in the SM, it is possible to use the measurement of Am^^ to obtain a second determination 
of fB^y/B~d and of fs^y/B^ (via ^s), thus reducing the error on the prediction for AMb^- 

The situation for Atub^ is different. From inspection of the standard fits of the unitarity 
triangle, it is clear that is always possible to choose p and rj such that the SM prediction agrees 
perfectly with the experimental central value. For this reason, in the numerical, analysis we 
just require the new physics contributions to Atub^ to be compatible with the experimental 
determination up to an uncertainty given by the lattice errors on f'B^y/Bd. A more correct 
analysis requires a simultaneous fit of the new physics contributions to Ek, Am^^, a^x^ and 
\Vub/Vcb\- See Ref. [75] for a general discussion of New Physics effects on Bg mixing. 

From the plots in Fig. |9] we see that Bq — Bq mixing data constraints are still much weaker 
than the corresponding constraint on Ek- 
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Figure 9: Plots a. m//± dependence of the T2HDM contributions to Am^^^^j in ps^ . Solid, 
dotted and dashed lines correspond to |^| = 1 and (tan Ph,'^^) = (30,0), (50,0), (50,7r/2), 
respectively. There is no appreciable dependence on C,'- The horizontal dashed lines are the 
experimental measurement. The blue band shows the theoretical uncertainties for the dashed 
line, similar bands can be drawn for the other curves. Plots b. Portions of the (tan/?//, mj^^i) 
plane excluded by Am^^^ ^j. The shaded area corresponds to ^ = 1. The dashed line show how 
this region changes for ^ = e*'^/^. Other choices of the phase yield in-between lines. 
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Figure 10: Plot a. mH± dependence of the T2HDM contributions to Aitib^^-^/ AniBf^^y See 



the caption in Fig. |9| Plot b. Excluded region in the ((p^,mH±) plane. The solid and dashed 
contours correspond to tan jSn =30 and 50, respectively. 



5.2.3 DD mixing 

The SM prediction for Sitld range between 10~^ ps~^ and 10~^ ps~^ and is completely dominated 
by long distance effects; in fact, the short-distance SM prediction has been calculated and 
reads [76,77] xd — 1-5 x 10~^ ps~^. The present experimental information on DD mixing 
parameters [78, 79], yields the following model independent determination of the DD mass 
difference [80]: Am^) = (14.5 ± 5.6)10~^ps~^. In the T2HDM very large effects are possible 
(of order 1% [8]), and there is the possibility that the actual D — D mass difference is entirely 
controlled by new physics short distance effects. In the numerics we require the new physics 
contribution to the Amo not to exceed the measurement. 
The D — D mass difference is given by 



Am 



D 



V uV* 

' uo ' cq 



VLL I r'VRR\ ir>VLL\ 
,D + '^l,D ) \V / 



+ 



SLL I ^SRR\ 
2,D "T ^2,D ) 



(57) 



where the matrix elements are 



LR 



{D\Q^^ 



\D) = —RmlflB!^^ 

\D) = -RmlflBl^^ 

\D) = -iRmlflB^^ 

\D) = iRmlflB^^ 



and R = {mD/{mc + mu)Y ■ In the numerical analysis we use B 



VLL 



B 



D 



(58) 

(59) 
(60) 
(61) 

(62) 

0.82 ± 0.01 [76] 
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Figure 11: Plot a. mH± dependence of the T2HDM contributions to Arrio- Solid, dashed 
and dotted lines correspond to |,^,^'| = 0.1, 0.2 and 0.5, respectively. We fix tsmPn = 50. The 
horizontal dashed line is the experimental upper limit. Plot b. Portion of the {C,C,',mH±) 
plane excluded by Am£,. The shaded area corresponds to tan/5/^ = 30. The dashed line to 
tan Ph = 50. 

(this value of the hat parameter B^, has been obtained from the lattice determination of 
i?£)(2GeV) [69,70]) and set all the other B parameters to 1. 

An approximate expression for the D — D mass difference is given by [8] 



Gtt' 



tan"^ [3h 



(63) 



The strong ^' dependence implies that, once B ^ tv data are imposed, no large deviations can 



be observed on Auid as can be seen from Fig 11 



5.3 B 



The branching ratio for the decay B — > rz/^ has been recently measured by the Belle [81] and 
Babar [82] collaborations 



B{B ^ TUr) = (l.79l°i^(stat);°i^(syst)) x 10"" [Belle] 
B{B Tv^) = (0.88+°:^?(stat) ± O.ll(syst)) x 10"^ [Babar] 

yielding the following world average 

BwAiB Tu^) = (1.31 ± 0.48) X 10"^ . 

The SM expectation reads: 



m. 



8tt 



1 - ^ fslVubl rB = (1.53 ± 0.38) x 10" 



m 



(64) 
(65) 

(66) 
(67) 



B , 
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tan(/3) '= 

Figure 12: Plot a. Portion of the (tan/3//, m^i) plane allowed by Rbtu for ^' = 0. Plot b. 
Excluded region in the {^',mH±) plane. The dotted, solid and dashed contours correspond to 
tanPn =10, 30 and 50, respectively. 



where we used the PDG world average \Vub\ 
surements only. The above result leads to 

Bwa{B 



R 



Btv 



(4.31 ± 0.3) X 10 ^ from direct tree level mea- 

> TUr) 



0.86 ±0.38. (68) 
= (3.68 ± 0.14) X 10"^ the prediction 



If we use the fitted value of the CKM angles {\Vuh 
reads Rbtv = 1-18 ± 0.50. The discrepancy between this determination of Rbtu and Eq. (68) 
is a manifestation of the conflict within the SM between the present determinations of Ki6 and 
sm{2p). 

In the T2HDM this process receives large tree level contributions via charged Higgs exchange: 

2 



tan^ Ph 



1 



m 



l-tan^/?//- ^ 



m 




(69) 



(70) 



The numerical impact of the constraint in Eq. (68) is very strong. In Fig. 12 d we show 
the impact of this constraint onto the (^',m//±) plane for various values of tan/?//. Since the 
experimental to SM ratio in Eq. (68) is smaller than 1, scenarios with ^' > are disfavored (see 
Eq. m). 



5.4 Time-dependent CP asymmetry in 5 ^ {J/i/j, (j),r]') Ks 

The time dependent CP asymmetry in the decay B f'Kg (/ = J/ip, 0, v') is given by 

(71) 



ajK - 


2 ImXfK 


l + \XfK\ 




(Mf/)* 


|Mf/| 
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where Fj^'* is obtained from Eq. (54) with obvious replacements. The effective Hamiltonian 
that controls the amplitude fK,) in the T2HDM is [50]: 



n 



V2 



.1=1 i=3 



(73) 



where 



Oi = 


{sl1,T''cl){cl1^T%l), 


(74) 


02 = 


isL7f^CL)icL'y^bL), 


(75) 


03 = 




(76) 


04 = 




(77) 


05 = 




(78) 


Oe = 




(79) 


OSQ = 




(80) 


0,Q = 




(81) 


0,Q = 




(82) 


0,Q = 




(83) 


Or = 


{cRhL) (slCr). 


(84) 



Tree level and one-loop charged Higgs diagrams contribute to the following matching conditions 
(we adopt the notation of Ref. [50]): 



R 



(^0) 



(/^o) 
(^0) 



(/^o) 



m 



iPRL)23 {Pj 



H \* 
RL) 22 



Vtb ml 



9slr 



DH{xth)+^CH{xth) (1 



K 



tan jSn 



3^2 



Vcb m]j^ 



tan' Ph 



tan' ISh 



(85) 
(86) 
(87) 
(88) 
(89) 
(90) 



where /xq ~ 0{mt), xth = m1/m'\j± and 
(^#l)33 (P#J32 tan^/?^ 



V* 



Vts \mtJ 



1-0.3^ 



V 30 ^ 



(91) 
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The functions Eh-, Dh and Ch can be found in Appendix A of Ref. [83]. Note that the results 
for the type-II 2HDM are recovered in the k — 1 hmit. 

Direct calculation of the anomalous dimensions involving the operator 0/j yields (in the 
notation of Ref. [50]): 

7^^ = -16, (92) 

7ff = (93) 

7{i°^ = 0(j^R,4). (94) 

(95) 

The large anomalous dimension of On implies a large impact of the running from the high-scale 
Hq ~ 0{mt) to the low-scale /i^ ~ 0(mb): 

CM = V-'-''''' Cnifio) , (96) 

where rj = Ois^iio) / ag^Hb) — 0.53. 

Let us first consider A{B^ J /ipKs). The impact of the QCD and electroweak penguin co- 
efficients is very small; hence the only T2HDM effect comes via the new operator Or. Adopting 
the naive factorization framework, the amplitude is proportional to: 

- 4 1 4 40 

A(SO^J/^ir,) (X -Ci(^fe) + -C2(/X,)+2C3(/X,)+20C5(^b) + 3C3Q(/i6) + yC5Q(/X6) 

-\Cn{^^,) , (97) 

where the Cr contribution receives a factor of —1/2 and 1/3 from Lorentz and color Fierzing, 
respectively. Note that the SM contribution (see for instance Ref. [84]) has been rewritten in 



the new operator basis Eq. (73). Direct calculation of this amplitude in the QCD-factorization 
approach [84] shows that the naive estimate is a fairly good approximation. 

The analysis of the amplitudes A{B^ — > (0, r]')Ks) is more complicated. In this case the QCD 
and electroweak penguin coefficients do play a leading role; moreover the magnetic penguin 
operator Og contributes as well. In the following we adopt the QCD-factorization analysis 
presented in Ref. [85]. After transforming the Wilson coefficients of Ref. [85] into our basis, we 
obtain the following expressions for the and t]' amplitudes: 

A{B° ri'Ks) oc (-0.04874 - 0.04905i) - (0.00658 - 0.00058i) 5C|"^(/io) 

-(0.00014 - 0.00002i) 5&^q\iiq) + (0.00711 - 0.00275z) 5C'i^\iiQ) 
+0.0007 Cj{iib) - 0.089 Cs{lib) + (0.03567 - 0.01087i) Cr{iiq) , (98) 

A{B'^^(j)Ks) oc (0.03262 0.00791i) (0.00963 - 0.00050i) (5C;5"\/Xo) 

+ (0.00044 - 0.00002z) 5d^Q\iiQ) - (0.00282 - 0.00013z) 5Cf °^(/io) 
-0.0004 C^{^^,) + 0.047 Cs{^b) - (0.01292 - 0.00092z) Cr{^io) . (99) 
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Note that the impact of the QCD and electroweak penguin matching conditions is suppressed 
by an order of magnitude with the respect to the leading contribution; in fact, the low-scale 
penguin Wilson coefficients are dominated by the tree-level coefficient C2 via the RGE running. 
Not surprisingly, the effect of the other tree-level operator {Or) on the running is also very 
large. 

The experimental measurements of these three asymmetries read [14]: 



a^K = 0.675 ±0.026, 
ar^'K = 0.61 ±0.07, 
a^K = 0.39 ±0.18. 



(100) 
(101) 
(102) 
(103) 



In Fig. [T3fi , we show the size of T2HDM contributions to the CP asymmetries in B ^ 
{J/ijj,r)',(j))Ks for some choice of input parameters. In Fig. 13 d, we show the portion of the 
(tan/?j:/, m^i) parameter space that is allowed by the present measurements of these asymme- 
tries. From the inspection of the figures we see that at the la level it is possible to reconcile the 
B ^ ip and B ^ t]' asymmetries in a quite wide region of the parameter space. The B ^ (p 
asymmetry, on the other hand, requires a too light charged Higgs. 



5.5 Ar,/r, 

The Bg-Bs width difference is given by 

AF, = -2 Tl^ cos(/3, ± Os) = -2{[Tl,]sM + STl^) cos(/3, ± 6^) 



12 



^(S,|Im|z J rf^xrKfr(x)Heff(0)}|5,) 



(104) 
(105) 



where, in the Standard Model, T^efr is the effective Hamiltonian that mediates the bottom 
quark decay and is dominated by tree level contributions. T2IIDM contributions affect both 
9s and F*2: the former has been already discussed in Sec. 5.2.2 the latter are induced by a 
new AB = 1 operator. The structure of the charged Higgs couplings in Eqs. (15) and (16), 



implies that the operator Or receives the largest contributions. In the numerics we will consider 
its effects together with the interference with the dominant Standard Model operator O2 (see 



Eq. ( 73 ) for the definition of the operators) 



Direct calculation of the T-product in Eq. (105) yields the following leading order expression 



for the new physics contribution to F 



SAT 



12 



^F'"'b,pole 

127r(2mBj 



1 - AZc Crr C2CRR 2 

2 4 2 c,pole b,pole 



(Q) 



(106) 
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Figure 13: Plot a. mH± dependence of the T2HDM contributions to Solid, dotted 

and dashed lines correspond to |,^| = 1 and (tan (3h, ^e) = (30, 0), (50, 0), (50, 7r/2), respectively. 
There is no appreciable dependence on ^' . The horizontal dashed lines are the experimental 
measurement. The blue band shows the theoretical uncertainties. Plot b. Portion of the 
(tan/3j:/, m//±) plane excluded by a(^^ri',(f>)K- The shaded area corresponds to ^ = 1. The dashed 
line show how this region changes for ^ = e*'^/^ (in the first plot, the region excluded is below 
the dashed line; in the second and third plots, it is above the uppermost dashed line and below 
the lowermost one). Other choices of the phase yield in-between contours. 



27 




Figure 14: Plot a. m^/i dependence of the T2HDM contributions to Ar^/Fs. We take |^| = 1 
and tan/?// = 50. Solid and dashed lines correspond to ip^ = and 7r/2, respectively. The blue 
band is the experimental 68% CL. allowed region. Plot b. dependence of the T2HDM 
contributions to AF^/F.,. We take |^| = 1 and tan/3// = 50. Solid and dashed lines correspond 
to m//± = 100 and 200, respectively. 



where the operators Q and Qs are 

Q = (bs)v-A(bs)v-A , 
Qs = (bs)s^p(bs)s-p 

and their matrix elements between 5, and states are 



iQs) 



m. 



r2 2 

[nib + nisY \ Nc 



2 - 



1 



with B = 0.87 ± 0.06 and Bg = 0.84 ± 0.05 [86]. After normalizing to the total B, 
obtain: 



AF. 



AF. 



SM 



+ tb, sat, cos{ps + Os) . 



(107) 
(108) 



(109) 
(110) 
width we 

(111) 



The SM prediction [87, 88] and the experimental result [25, 89] read: 

■AF. 



F. 
AF. 



SM 



exp 



0.147 ±0.060 
0.27 ±0.08 . 



[113) 



28 



5.6 CP asymmetry in flavor speciflc B decays 

The CP asymmetry in flavor specific Bg decays (only the decays Bq f and Bq ^ J are 
allowed) is given by: 



A 



(<?) 

SL 



T{Bq{t) ^ /) - V{Bq{t) f) 
T{Bq{t)^ f)+T{Bq{t)^f) 



Im- 



12 



M 



(<?) 

12 



(114) 



From the discussion in Sec. 



5.5 



it follows that the T2HDM effects on r^2^ are negligible, hence we 
will consider only box diagram contributions to M^f ■ Adopting the standard parametrization, 



mS/m[%,, = r2exp(2^^,) we get: 



A 



SL 



Im 



12 



COS 29 „ 



-Re 



12 



sin 26n 



12 / SM 



[115) 



where [86,90] 



(r(?/M(? 



(rg/M: 



SM 



SM 



-(40±16)-i(5±l) 



X 10 



-4 



(40± 16) +i (0.22 ±0.04) 



X 10" 



(116) 
(117) 



In the notation of Sec. 
the SM predictions: 



5.2.2 



we have exp(2z6'g) = F^^'^ / So{xt) ■ From Eqs. (116 117) we read 



SM 
SM 



(5 ± 1) X 10"^ 
(0.22 ±0.04) X 10"^ 



:ii8) 

:il9) 



Unfortunately the experimental errors on these asymmetries are at least an order of magnitude 
larger than the SM expectations [91-96]: 



[A'A] 



exp 



exp 



;il ±55) X 10"^ 
;-80± 110) X 10"^ . 



(120) 
(121) 



In Fig. 15 we show the size of the T2IIDM contributions to these asymmetries for large tan (5h- 



5.7 Neutron EDM 

The effective Hamiltonian that encodes charged Higgs contributions to the neutron EDM is 

e 



cflf 



q=u,d 



167r^ 



(122) 
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400 



Figure 15: tp^ and m//± dependence of ^5^' fo'^ 1^1 = 1 ^'^^ tan (3h = 50. The dashed hnes 
are the la SM expectation. The blue band is the theoretical error. 
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Figure 16: Constraint that the neutron EDM puts on the T2HDM parameter space. The 
shaded area is excluded at 90% CL. The solid and dashed lines correspond to mH± = (200, 1000) 
GeV. 
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where 



^ ^ I (PlrUPrl):! ^ (!I!L\+...{ (123) 



and the function 5 is given in Ref. [36]. The dots stand for terms that do not contribute to 
the imaginary part of the Wilson coefficients. In the chiral quark model, the neutron EDM is 
related the valence quark EDM's, dy, and dd, via [97]: 

dn = \{^dd- dy) T]^ , (125) 
= yI^MC,), (126) 

where r]^ ^ 1.53 is the QCD correction factor. Approximate formulae for the up and down 
quark contributions to the neutron EDM are (in units of e cm): 



^rj^dd = 10"''^(^^)'(;^j [-(6.5 ±0.3) cos ^^-(15.9 ±0.7) sin (127) 

where the uncertainties come from varying in the (200 — 1000) GeV range. Taking 

into account that the 90% CL. experimental upper bound on the neutron EDM is [25] 6.3 x 
10~^^ e cm, it is clear that the T2HDM parameter space is constrained only for ^' ^ 0. Note 
that the huge enhancement in du comes from S13 oc mctan Ph^' /mw that is not suppressed 
either by Vub or, for large tan/5//, by the charm Yukawa. 

In Fig. [16] we show the impact of the present upper bound on the T2HDM parameter space. 



2 



5.8 CP asymmetries in B ^ K tt and ^ K 7r+ 

The direct CP asymmetries in the decays K^n^ and — > R-tt^ can be calculated 

(albeit with large errors) in the QCD factorization approach [31]: 

Acp{B-^K-n') = (7.11}:™^:^:?)% (129) 

Acp{B'^K-n^) = (4.5^}:i™|iL^) % (130) 

where the first error corresponds to variation of the CKM parameters, the second and third 
errors refer to uncertainties in the hadronic parameters used in the calculation and the fourth 
error refiects additional uncertainties caused by the breakdown of the factorization ansatz (that 
result in endpoint singularities regulated in terms of two extra complex parameters). Because 
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150 200 250 300 350 400 450 500 -150 -100 -50 50 100 150 



Figure 17: mH± and ip^ dependence of the T2HDM contributions to AAcp. Solid, dashed 
and dotted hues correspond to (tan/^^, = (50,1), (35,1) and (50,2), respectively. The blue 
band is the experimental 68% CL. allowed region. In plot (a) and (b) we fix ip^ = —50° and 
mH± = 200 GeV, respectively. 



of a high degree of correlation, most of these errors cancel when considering the difference 
between these two asymmetries. From Ref. [31], we see that this difference lies in the range 
[0.5,3.3]. Using the formulae presented in Ref. [98] and updated numerical inputs, we find 

AAcp = Acp{B- ^ K-Ti'^) - Acp{B° ^ K-7T+) = 2A{1±0.5) , (131) 

where the 50% error refiects the uncertainties studied in Ref. [31] and the possibility unusually 
large power corrections (scenarios S1-S4 of Ref. [31]). This SM estimate has to be compared 
to the present experimental results [14]: 

Acp{B- ^ R-TT^) = (4.7 ±2.6)% (132) 
Acp(5° ^ ir-7r+) = (-9.7 ±1.2)% (133) 
AAcp = (14.4 ± 2.9) % . (134) 

The effective Hamiltonian responsible for T2HDM contributions to AAcp and the relevant 



matching conditions have been given in Sec. 5.4 For completeness we point out that other 



approaches to the calculation of AAcp (see, for instance, Ref. [99]) are consistent with the 
QCD-factorization results. Also note that it might be possible to accommodate the present 
experimental results in models in which the color-suppressed tree contribution is unusually 
enhanced [100]. 

In Fig. [17] we show the size of T2IIDM contributions to AAcp- Unfortunately, for m//± 
larger than 300 GeV we do not find any sizable effect. 
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6 Observables: the neutral Higgs sector 



6.1 (^-2V 

The anomalous magnetic moment of the muon receive potentially large contributions at the 
2-loop level via the Barr-Zee mechanism [101, 102]. These diagrams are able to account for the 
large discrepancy between the SM prediction and the experimental measurement only for very 
light pseudo-scalar mass {uia < 100 GeV) [103]. 

6.2 Ap 

T2HDM contributions to Ap depend on both the neutral and charged Higgs bosons. For any 
choice of the charged Higgs mass it is possible to find region of the (m^, uih and an) parameter 
space for which the corrections to the p-parameter are in agreement with the experimental 
bounds. 

6.3 Z ^bb 

Charged and neutral Higgs contributions to the effective Z — b — b coupling affect both the ratio 
i?6 = r{Z — > bb)/T{Z hadrons) and the forward-backward asymmetry A^. The present 
experimental results are [25]: 

Rb = 0.21629 ± 0.00066 (135) 
Ab = 0.901 ± 0.013 , (136) 

where the value for Ai, has been obtained by combining the direct and indirect measurements 
(from Aft = 4/3 A^^^/A^). The SM fits for these two observables read: Rf^ = 0.21586 and 
A^^ — 0.935. In particular note that Af, shows a deviation of about 2.5(7 from the SM predic- 
tion. 

The T2HDM contribution to the effective Z ^ bb vertex can be easily extracted from the 
results of Ref. [104]. These corrections depend on both the charged and neutral Higgs sector of 
the T2HDM. Unfortunately we do not find any sizable effect in the portion of the parameter 
space allowed by the other constraints. 

7 Summary & Outlook 

Thanks to the spectacular performance of the two asymmetric B-factories, intensive studies in 
the last few years have demonstrated that the CP and flavor violation observed in B and K 
physics is described by the standard CKM mechanism to an accuracy of about 15%. A very 
interesting result, potentially one of the most important discovery made at the B-factories, is 
that the time dependent CP-asymmetries in penguin dominated modes do not seem to agree 
with the SM expectations. At the moment these deviations are in the 2.5-3.5(7 ranges. Since 
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these modes are short-distance dominated, they are very sensitive to presence of beyond the 
Standard Model phases. Taking seriously this pieces of data is suggestive of sizable contri- 
butions from a non-standard phase. For the sake of completeness we also mentioned several 
other measurements that display a significant deviation from the SM, such as difference in the 
CP asymmetry between K^ir' and K^i:^, the (g-2) of the muon and the forward-backward 
asymmetry m Z ^ hh measured at LEP 

As an illustration of a new physics scenario that may account for the observed deviations from 
the SM we have presented a detailed study of the two Higgs doublet model for the top quark. 
The model is a simple extension of the SM which in a natural way accounts for the very heavy 
top mass. We view it as an interesting low energy effective model that encompasses some of the 
important features of an underlying framework of dynamical electroweak symmetry breaking. 
Of course, the deviations seen in B decays and other flavor physics, may also be accountable 
by many other extensions of the SM; for example supersymmetry [10], a fourth family [11], a 
Z-penguin [12] or warped flavor-dynamics [13]. Obviously, the main features of any extension 
of the Standard Model that is to account for the experimental deviations in B-physics and 
other flavor physics that are discussed here are that there have to be new particles in the ~ 
300 GeV to ~ few TeV range and associated with these one needs at least one new CP-odd 
phase. Distinguishing between various scenarios or nailing down the precise structure of some 
other models responsible for the deviations that we have discussed will undoubtedly require 
much more experimental information. In any case if the hints from the B-factories are really 
true, we will certainly witness a truly exciting era in Particle Physics as the LHC starts its 
long awaited operations in 2008. It should also be abundantly clear that infusion of precise 
information from low energy flavor measurements will be crucial in interpreting the findings at 
the terascale energy. 
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